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A simple experimentally accessible realization of current rectification by molecules (molecular films) 
bridging metal electrodes is described. It is based on the spatial asymmetry of the molecule and 
requires only one resonant conducting molecular level (n orbital). The rectification, which is due to 
asymmetric coupling of the level to the electrodes by tunnel barriers, is largely independent of the 
work function difference between the two electrodes. Results of extensive numerical studies of the 
family of suggested molecular rectifiers HS-(CH2)m-C6H4-(CH2) n -SH are presented. The highest 
rectification ratio ~ 500 is achieved at m = 2 and n = 10. 

PACS numbers: 85.65.+h 



I. INTRODUCTION 

The concept of a molecular rectifier started the field of 
molecular electronics (moletronics) almost thirty years 
ago [ID. At present, moletronics is an explosively grow- 
ing field of experimental and theoretical activity. The 
molecular rectifier (MR) is still one of the central objects 
of this research. Even the simplest future application, 
such as molecular memory, will require high quality MR 
with sharp voltage thresholds, large current rectification 
ratios, small time constants, large breakdown voltages 
and so on. MR demonstrated so far ]^-[4| constitute an 
impressive proof of principle but their electrical parame- 
ters are in fact very poor. The situation prompts further 
investigation of molecular rectification. 

Until recently, theoretical analysis of MR had been 
largely limited to the Aviram and Ratner donor- 
insulator-acceptor (D — a — A) mechanism jjj. In this 
design, the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO) are 
confined to two different parts of the rectifier, D and A re- 
spectively. The insulating bridge er prevents the orbitals 
from "spilling off" to the other part. If such a molecule is 
placed between two metal electrodes, the current-voltage 
characteristics of the junction is expected to be highly 
asymmetric. At a particular voltage applied in the posi- 
tive direction, the Fermi level of the electrode on A side 
aligns with the LUMO, and the on D side aligns with 
the HOMO. At this voltage, the current rises sharply be- 
cause the electrons can now be loaded on the LUMO, 
then tunnel inelastically through the a to the HOMO 
and the escape into the second electrode. In the opposite 
direction, a similar process does not occur until a much 
higher applied voltage. 

A different mechanism of molecular rectification was 
described by Ellenbogen and Love ||. It is based on an 
energy mismatch between two conducting levels localized 
on different parts of the molecule. Within the D — a — A 
framework, such levels could, for instance, be the LU- 



MOs of D and A. Under external bias, the levels shift 
because of the electric field. In the forward direction, the 
conducting levels align at some voltage, and facilitate 
resonant transport of electrons between the electrodes. 
In the reverse direction, the levels move away from each 
other and the current remains small. Hence molecular 
rectification. 

We note that both described mechanisms require two 
electroactive molecular levels and a fine balance between 
potential drops inside the molecule and on the molecule- 
electrode interfaces. To the best of our knowledge, nei- 
ther of the mechanism has been realized experimentally. 

Recently, an interesting case of molecular rectifica- 
tion was discussed by Krzeminski et al. They 
studied theoretically current rectification in Langmuir- 
Blodgctt (LB) films of 7-hcxadccylquinolinium tri- 
cyanoquinodimethanide (C16H33Q-3CNQ), which was 
previously studied experimentally |^|,^|. That molecule 
was initially considered to be a possible implementation 
of the Aviram-Ratner mechanism because its active part 
Q-3CNQ comprised a donor and an acceptor. (The long 
insulating tail C16H33 was added to help form good LB 
films.) It was later realized J5J that the Q-3CNQ is, in 
fact, a D — 7r — A molecule and is unlikely to implement 
the Aviram-Ratner mechanism. The analysis of Ref. |J 
confirmed that the 7r bridge does not sufficiently isolate 
D and A to keep molecular orbitals localized on cither 
D or A. Instead, the orbitals were delocalized over the 
entire Q-3CNQ unit. Krzeminski et al. have attributed 
the observed rectification to asymmetric position of the 
LUMO and HOMO with respect to the Fermi levels of 
the metal and to "asymmetric profile of electrostatic po- 
tential across the system" . 

The purpose of this paper is to point out that there 
is, indeed, a simple and general mechanism of molecu- 
lar rectification where a single electroactive unit is posi- 
tioned asymmetrically with respect to electrodes and the 
HOMO and LUMO are positioned asymmetrically with 
respect to the Fermi level. However, this mechanism does 
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not require such complex molecules as Ci6H 33 Q-3CNQ. 
What is needed is just one conducting molecular level 
placed closer to one electrode than to the other. Since 
most of the applied voltage drops on the longer insulat- 
ing barrier, the conditions for resonant tunneling through 
the level are achieved at very different voltages for the 
two opposite polarities. By changing the lengths of the 
insulating barriers, the rectification ratio can be system- 
atically changed. This mechanism can be realized by 
relatively simple molecules. For instance, the conduc- 
tive level could be supplied by a benzene ring. Also, the 
molecules can be shorter, produced by e.g. self-assembly 
instead of LB deposition, and therefore more conductive. 
Indeed, the currents reported in |3) were of order 10~ 17 
A/molecule, which almost rules out any practical appli- 
cation of such MR. This is obviously a result of having 
the long aliphatic tail C16H33. There is an experimen- 
tal evidence that even much shorter alkane chains, like 
C12, are very resistive and transport there proceeds by 
hopping processes rather than tunneling ||. We show 
below that With simpler conducting and shorter insulat- 
ing units, MR can achieve rectification ratios in excess of 
a hundred while remaining fairly conductive. 

In the foregoing sections, we discuss the present mech- 
anism of molecular rectification in more detail. We have 
performed a numerical analysis of the mechanism by cal- 
culating the current-voltage (I-V) characteristics of the 
family of prototype molecular diodes HS-(CH2) m -CeH4- 
(CH2) n -SH. We have found that rectification ratios of 
> 100 are achievable with such a design. 

We mention for completeness that any asymmetric 
electrode-molecule-electrode junction should in princi- 
ple produce asymmetric currents at large enough volt- 
ages. Such a high- voltage asymmetry due to unequal cou- 
pling to the electrodes [y,[10| or to an asymmetric central 
molecular unit jnj have been observed experimentally 
and discussed theoretically. However, the current rectifi- 
cation ratios in these studies were of order unity, which 
is clearly insufficient for practical applications. Zhou et 
al reported rectification in a molecular monolayer of 4- 
thioacetylbiphenyl m. In that experiment, the two elec- 
trodes were different, with different work functions, and 
different connections to the molecules. The experimental 
data were interpreted as the standard thermionic emis- 
sion through an asymmetric barrier. 



II. A MOLECULAR DESIGN FOR CURRENT 
RECTIFICATION 

Our molecular rectifiers with asymmetric tunneling 
barriers consist of five structural parts, see Fig. |. The 
end groups provide contact to the two electrodes. For 
instance, these may be thiols that could self-assemble on 
the surface of gold or silver, carboxyl groups in case the 
Langmuir-Blodgett technique is used, and so on. The end 
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FIG. 1. Top: Schematic structure of an asymmetric tun- 
neling barrier molecular rectifier. Bottom: The (2,6) member 
of the MR family HS-(CH2) m -C 6 H4-(CH 2 )„-SH. 



groups may even be absent altogether if the application 
does not require them to provide better contacts. The 
three inner parts of MR are the central conjugated group 
C and two insulating barriers 1l and 1r. The purpose of 
such a construction is to provide an electronic level local- 
ized on C, with energy not very different from the Fermi 
energies of the electrodes. In most cases this level will 
be the lowest unoccupied molecular orbital (LUMO) of 
the molecule. Then the electron transmission probability 
should be resonant near the energy of the LUMO. 

The energy diagram of MR is shown in Fig. [| The two 
main parameters that determine the rectification proper- 
ties of MR are the energy difference between the LUMO 
and the Fermi energy of the electrodes A = W — A, and 
the ratio of the voltage drops on the right and left insu- 
lating parts 77. Under the assumption that the polariz- 
ability of C is much larger than that of the Is, the voltage 
drops on the barriers are proportional to their respective 
lengths L r ight and Li e ft, and 77 w Lright/^left- For simplic- 
ity, we assume that the two electrodes are similar, or have 
similar work functions. Therefore there is no contact po- 
tential difference and the electric field on the molecule is 
zero at zero applied bias. The minor complications that 
may arise from the contact potential difference will be 
briefly discussed in Section [v|. Let us consider the op- 
erating principle of MR. We shall follow the convention 
that the right electrode is always grounded and the right 
Fermi energy can be used as the reference energy. Then 
a positive potential applied to the left electrode shifts its 
electronic levels to lower absolute energies. Thus under 
positive bias, the Fermi energy of the left electrode goes 
down. Due to the non-zero electric field, the LUMO will 
be dragged down too. The energy shift of the LUMO 
A-Elumo is determined by the parameter r\. Under our 
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FIG. 2. The basic principle of rectification by asymmetric 
tunnel barriers. W is the work function of the metal, A is the 
electron affinity of the molecule, A = W — A. (b) A molecule 
with A <W and with different lengths of the insulating bar- 
riers is connected to two metallic leads, (c) Under forward 
bias, the current rises when the right Fermi level aligns with 
the conducting molecular level, (d) Under reverse bias, the 
current rises when the left Fermi level aligns with the con- 
ducting molecular level. Since most of the total voltage drops 
on the right insulation barrier, Vr > Vf- 



assumptions, the field across C is very small, and the po- 
tential drops on the two barriers, £/i e ft and [/ r i gri t have to 
sum up to the total external bias V, f/i e ft + [/right = V. 
Then it is easy to show that 
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where q is the elementary charge. A sharp increase in 
current is expected when the LUMO lines up with the 
right Fermi level, that is at a£lumo = A. The corre- 
sponding forward voltage follows from Eq. (pi) 
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Under reverse bias, both the left Fermi level and the 
LUMO go up in energy. The current will turn on when 



the left Fermi level lines up with the LUMO, when 
A = [/i ft = V — C/ r ight- This corresponds to a reverse 
voltage 
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Clearly, the forward and reverse voltages are different, 
their ratio being Vr/Vf — V- At a large 77, the two 
thresholds differ, and there is a voltage window Vf < 
\V\ < Vr within which there is a substantial current in 
the forward direction and almost no current in the reverse 
direction. Hence, a strong rectifying effect is expected. 

Different parts of MR affect different properties of the 
I-V characteristic. Parameter A determines the overall 
scale of the forward and reverse voltages. A itself could 
be systematically changed by using different materials 
for the electrodes and different conductive parts C. The 
work functions of the electrodes used in moletronics stud- 
ies vary from 5.6 eV (Pt) to 4.0 eV (Si), enabling almost 
continuous adjustment of A. (Obviously, any changes 
of the electrode material have to be compensated by re- 
spective changes of the side groups of MR.) The use of 
different materials for the two electrodes provides fur- 
ther tuning of A. Indeed, the contact potential changes 
energy of the LUMO by some additional amount. [Cf. 
Eq. (Q) and use the contact potential in place of bias V.] 
The conclusion from this qualitative analysis is that A is 
a parameter under the control of a designer of molecular 
rectifiers. What absolute value of A, large or small, is 
optimal from the electronic viewpoint is not clear a pri- 
ori. Increasing A widens the rectification window and 
improves the stability of the device. At the same time, 
larger As imply higher operating voltages, power dissi- 
pation and other unwelcome consequences. 

The shortest of the two insulating barriers controls the 
width of the conducting level and consequently the sharp- 
ness of the current increase in the vicinity of Vp and 
Vr. Indeed, the transmission probability T(E) through 
a molecular level, which defines the conductance of the 
molecule G(E) oc T(E) according to standard Landauer 
formalism [14| , is well described by the Breit-Wigner for- 
mula 
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where Emo is the energy of the molecular orbital avail- 
able for resonant transmission (LUMO in the present 
case), and T^r are the partial widths associated with 
coupling to the left and right electrodes, respectively, 
see analysis in |lq| . The "golden rule" estimate gives 
Tl.r ~ t\ r/Dl,r, tL,R being the effective coupling ma- 
trix clement and Dlr the electron bandwidth in the elec- 
trodes. We see that a steep rise in current occurs when 
one of the electrode Fermi levels rises to line up with 
EmOi which is exactly what is required for large rectifi- 
cation ratios. If the C component of the molecule is too 
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close to a metal, the conducting level is so broad that the 
transmission probability will be substantial at all ener- 
gies. Indeed, I\ oc e~ 2KLi , where k oc (-Efcamer — Ep) 1 ! 2 is 
the tunneling attenuation coefficient. Here -Ebarricr is the 
energy of of the conducting level of the insulating bar- 
rier, which lies higher that the LUMO of the conducting 
unit C. For alkane chains on gold, (-Ebarrier — Ep) ks 4.8 
eV H . Thus the width of the level is exponentially sen- 
sitive to the distance from the electrode. If the level is 
substantially broadened, there will be little difference be- 
tween the currents in the forward and reverse directions. 
Sufficient insulation of the conducting molecular orbital 
from both electrodes is an essential feature of the present 
rectification mechanism. 

Finally, the longest of the two barriers controls the 
anisotropy of the I-V characteristic and the overall am- 
plitude of the current. As was discussed above, Vr/Vf ~ 
-bright /Li e {t ■ Therefore the longer the second barrier the 
larger the Vr/Vf ratio and the better the diode. At the 
same time, the current goes down exponentially with the 
barrier length. Indeed, the current is 
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when the resonance falls into the "window" between the 
lowest and the highest Fermi levels in the leads. The 
current falls off exponentially with the thickness of the 
thicker of the insulating barriers L > . Therefore, we are 
facing a typical trade-off problem but in its worst form. 
The asymmetry of MR improves linearly but the resis- 
tance grows exponentially with the length of the bar- 
rier. Nonetheless, we found that substantial rectification 
is achieved at the current levels that make electronic ap- 
plications of MR practical. This is detailed in the next 
sections. 



III. CALCULATION DETAILS 

In order to test the proposed mechanism of molecular 
rectification, we calculated the I-V characteristics of a 
prototype family of molecular diodes HS^CR^n-CeHzi- 
(CH 2 ) m -SH sandwiched between two gold electrodes. 
The five functional parts of MR are constructed from 
the standard chemical groups used in moletronic stud- 
ies. The end thiols -SH chemically attach the molecules 
to gold, the middle benzene ring provides a conducting 
level (LUMO) at energy -Elumo = —3.5 eV (with respect 
to vacuum) , and the insulating barriers are made of sat- 
urated hydrocarbon units -CH2-. Individual members of 
the family are parameterized by the two integers (n, m) 
which are the respective numbers of the -CH2- groups on 
both sides of the molecules. 
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FIG. 3. A monolayer of molecular diodes between two 
semi-infinite fee electrodes. Only one layer of electrode atoms 
is shown. The solid line encloses the atoms included in one 
molecular complex. 



The electrode-MR-electrode junction is shown in 
Fig. 0. The equilibrium structure of isolated molecules 
was obtained through total energy minimization with the 
density- functional program Spartan | |l3fl . An important 
factor affecting the electrical properties of the junction 
is the local geometry of the molecule-electrode contact. 
It is generally accepted that, upon self-assembly on gold, 
thiols lose the end hydrogens and then bind directly to 
the gold atoms. The preferential binding position of sul- 
fur atoms on the Au(lll) surface is at the apex of the 
pyramid with a triangular base of gold atoms (hollow po- 
sition). Also possible is the top position, where the sulfur 
is positioned directly above one of the gold atoms |l^] . In 
this paper, we assume the hollow binding position of the 
sulfur. In the hollow position, the end S states strongly 
hybridize with states on gold atoms, and the resulting 
current is less dependent on exact molecule-contact ge- 
ometry in comparison to the top position jOI. Start- 
ing with the equilibrium molecular structure we remove 
the end hydrogens and replace them with two clusters 
of three gold atoms. The clusters form cr-bonds to sul- 
furs instead of the lost hydrogens. The primary molecule 
axis is roughly perpendicular to the planes of the gold 
triangles. The resulting molecular complexes (three gold 
atoms - dithiolate - three gold atoms) are organized in 
a periodic two-dimensional film which is placed between 
two (111) surfaces of semi-infinite gold electrodes. The 
film is commensurate with the Au(lll) surface but its 
primitive surface cell is four times larger. 

In order to obtain the I-V characteristic of the device 
we use the Landauer formulation of quantum transport 
|l4| together with a semi-empirical tight-binding param- 
eterization of the Hamiltonian matrices for the molecule 
and the electrodes. The off-diagonal matrix elements of 
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TABLE I. Tight-binding parameters of atoms used in 
transport calculations. 



the molecule are taken from Harrison 
Slater-Koster linear combinations [jl6| 
elements W a py 



fL5J. They all are 
of the four basic 



W afh (i - i) = 7 -62^1 eV, 



(6) 



where tiy is the distance between atoms 



i and j, and 

Vssa = -1.40, ry spo - = 1.84 r\ ppa = 3.24, and r] pp7r = 
—0.81. The diagonal matrix elements, that is on-site en- 
ergies of core atomic orbitals, are also adopted from 
|l5| , while the on-site energies of valence orbitals are cal- 
culated from atomic affinities A; and ionization poten- 
tials L as discussed in H 



A, 



Si = — 
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The affinity and ionization potential determine the 
atomic Hubbard parameter Ui as well: 



(8) 



Table | lists the molecular tight-binding parameters used 
in this paper. 

Our parameterization for the electrodes is based on the 
parameter sets of Papaconstantopoulos Although 
our calculational procedure enables us to treat the leads 
with the full set of s, p, and d orbitals JT^j2TJ], for the 
purposes of this paper the details of the electrode elec- 
tronic structure are not important. We therefore choose 
to work with the s component of the band structure only. 
The W ssr7 parameter of Ref. JTJ for gold is -0.909 eV. 
We found, however, that this value results in too nar- 
row a band with the high-energy edge being not very far 
from the conducting LUMO of MR. In order to separate 
molecular rectification from the specific effects associated 
with gaps in the leads' density of states, we also used the 
double value, W ssa — —1.818 eV. The same value of W ssa 
is used for matrix elements inside the electrode and be- 
tween the electrode surface and atomic triangles that are 
part of the molecules. A comparison of the results ob- 
tained using the two values of W ssa will be given in Sec- 
tion [f^. The second important parameter is the Fermi 
energy of the electrodes (that is minus the work function 
of the lead material) . The effectiveness of our rectifica- 
tion mechanism depends critically on Ep. Therefore we 



choose to study a series of Ep which imitates the effect 
of using different electrodes. In particular, we studied 
the six values E F = -4.0, -4.2, -4.5, -5.0, -5.2, and 
-5.5 eV. 

Calculation of current is performed utilizing a multi- 
step procedure. First, an isolated molecular complex 
(three Au atoms - dithiolate - three Au atoms) is treated 
as follows. The molecular Hamiltonian -ff mo i is con- 
structed from the parameters described above. Then the 
molecular wave functions and their energies are found by 
diagonalizing the matrix E — H mo \. From the wave func- 
tions, the average number of electrons on each atom is 
calculated: 



(9) 



where the index n numbers the molecular orbitals and on 
the atomic orbitals that belong to atom i. After that the 
diagonal matrix elements of H ma \ (the onsite energies St) 
are recalculated as 



£■ = £i + Ui(qi - Zi 



(10) 



where Zi is the atomic charge of the i-th atom, which in 
turn changes the charges qi on the atomic sites. Then 
the procedure is repeated until the charges qi converge. 
The converged charges define the final position of the 
molecular levels and the molecular Hamiltonian to be 
used in the transport calculations. 

In the second step, the semi-infinite electrodes are 
solved. In our approach, there are two global quantum 
numbers, total energy E and momentum parallel to the 
surface k|| . Fixing k|| converts the semi-infinite three- 
dimensional problem into a semi-infinite one-dimensional 
problem with a finite dimensional basis. Notice, however, 
that the Hamiltonian of the one-dimensional wire explic- 
itly depends on k|| . Thus transport problems at different 
ky are not equivalent. Moreover, the dependence on this 
quantum number could be substantial. It is important 
therefore to consider a grid of k|| and the grid should be 
as dense as possible. Then we follow the procedure of 
Ref. PQ| , modified for the presence of an arbitrary ori- 
ented surface. For each E and k|| , we solve the channel 
problem, i.e. find all the Bloch vectors with both real 
and complex k z vectors. The former correspond to open 
conducting channels while the latter to evanescent chan- 
nels. The wave functions and k z values of the channels 
are used to construct the surface Green's functions of the 
electrodes. Note that no energy or momentum integra- 
tion is required with this method. 

On the next step, the molecule is eliminated from the 
picture by first solving the Schrodinger equation for the 
molecular wave function and then substituting the solu- 
tion into the Schrodinger equation for the surface atoms 
of the leads. As a result of such a procedure, a matrix 
operator V appears that directly couples states on the 
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left wire with the states on the right wire. V is essen- 
tially the inverse of E — H mo i convoluted with the matri- 
ces that describe the molecule-electrode coupling. There- 
fore, it contains all the information about the molecule. 
In particular, V has poles at the energies of the molec- 
ular orbitals. Knowing V and the Green's functions of 
the free electrodes, the full electrode Green's function is 
found from the Dyson equation [EfJ 



g = (Go- 1 - vy 



(11) 



where Go is the block-diagonal matrix in which the upper 
left corner is the free surface Green's function of the left 
wire and the bottom right corner is that of the right wire. 
Note that the size of all the matrices involved in Eq. ( pd| ) 
is equal to the combined dimension of the left and right 
electrode surfaces. In the present study, each surface unit 
cell has four states, therefore all the matrices are (8x8). 

Since the Green's function solves Schrodinger's equa- 
tion for non-coincident spatial arguments, the transmis- 
sion coefficients t nn '(E,\s.\\) between the open channels 
can be found by multiplying G by certain projector vec- 
tors P. Each projection recovers an incident wave in 
one of the open channels in a particular electrode. The 
components of the vector GP determine the transmission 
and reflection amplitudes in all the open channels of both 
electrodes 20 . The contribution to current is found by 



squaring the transmission coefficients and summing over 
the open channels of the receiving wire. The total cur- 
rent is obtained by summing over all the incident chan- 
nels available, then over k|| , and integrating over energy 
E: 



In r E F +qV/2 



lEp-qV/2 



T{E)=Y J Vnn>(E^)\ l . 



(12) 

(13) 



In the present study, we have used a grid of 64 k|| points 
uniformly distributed over the two dimensional surface 
Brillouin zone. 

An important and technically difficult problem is to ac- 
count for the additional charge that the molecule acquires 
due to the mismatch of its equilibrium chemical potential 
with the Fermi energy of the electrodes. All the molec- 
ular components of the scattering wave functions could 
be found from the same Green's function (^). Squaring 
them and integrating over energy, one finds new charges 
on the molecule, much like Eq. (|^). Then the new on-site 
molecular energies are computed according to Eq. (10). 



Thus in principle, the entire calculational procedure de- 
scribed above (except the final calculation of current) 
has to be repeated many times until the charges con- 
verge. Moreover, this self-consistent procedure has to 
be performed anew for every value of the external bias 
voltage V, i.e. out of equilibrium. The charges in this 
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FIG. 4. The bulk electrode density of states for the two 
values of W asa . The data are collected from a (8 x 8) grid of 
k|| points, the one-dimensional square-root singularities being 
visible. The vertical solid line indicates the position of the 
Fermi level. Inset: the same density of states but computed 
on a (40 x 40) grid of ky points. 



non-equilibrium case can be calculated with the use of 
the procedure described in |i~9fl . 

The additional charge transfer is most important for 
strong coupling between the molecule and the wires. In 
this case, the molecular levels broaden significantly so 
that their tails can accommodate significant additional 
charge. Luckily, our rectification mechanism works best 
for narrow molecular resonances. As long as the Fermi 
energy remains below the conducting molecular level, the 
additional charge is negligible. For this reason, we do not 
perform the self-consistent calculation of charge in this 
paper. What we do take into account is the linear shift 
of the on-site atomic energies due to the external electric 
field due to the bias voltage. 



IV. RESULTS 

We begin by presenting our numerical results for the 
electrode density of states N(E) (DOS), see Fig. ||. 
Clearly visible are the square root singularities, which are 
characteristic of one-dimensional conductors. Recall that 
fixing ky renders the wire one-dimensional and the three- 
dimensionality is restored upon summation over the in- 
finite number of kii points. The comparison of the main 
panel and inset in Fig. |i| indicates that this is indeed the 
case. Our choice of 64 kii points is a reasonable com- 
promise between the accuracy and time for the transport 
calculations. The DOS affects the I-V characteristic of 
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FIG. 5. Transmission probability through molecular 
diodes Au3-S-(CH2)„-C 6 H4-(CH2)6-S-Au3, n = 0,1,2. The 
peak corresponds to transmission through the LUMO located 
on the benzene unit. 

the diode in two ways. First, it changes the effective 
resonance width since T oc N(E). Thus, in the case of 
W ssa — —1.818 eV, the resonances are expected to be 
twice as narrow as for W SS(J — —0.909 eV (N oc 
Secondly, gaps/edges in the density of states could block 
clastic tunneling completely, resulting in a strong nega- 
tive differential resistance effect. This possibility is dis- 
cussed in Subsection (iv) below. 

We investigated four major effects of the spatial and 
electronic structure of the junction on its I-V character- 
istic. 

(i) Effect of the width of the molecular level on the 
rectification property of the junction. Within our mecha- 
nism, good rectification requires narrow molecular reso- 
nances. If the levels are narrow, the current rises steeply 
upon reaching the threshold voltage. In other words, 
there is very little current before the resonance is reached 
but some finite current after, which can produce a high 
rectification ratio. In the case of a broad resonance, the 
Lorentzian tails of the transmission function provide a 
sizeable current even before the nominal threshold volt- 
age is reached. The rectification ratio should be small in 
this case. The width of the resonance is controlled by the 
thickness of the shorter insulating barrier. In Fig. ^ we 
compare the transmission function for the three molec- 
ular rectifiers which have no, one, and two insulating 
-CH2- groups on the short side. One can see that in the 
first two cases the resonance has a sizeable width of 0.2 
eV (full width at half maximum), which implies strong 
coupling to the electrode. This is more or less obvious 
for n — 0, when the sulfur atom is directly attached to 



the benzene. The lone electron pair of the sulfur overlaps 
with the 7r electrons of the ring, resulting in a molecu- 
lar orbital distributed almost evenly over the ring and 
the sulfur. Since the sulfur is directly coupled to the 
electrode (albeit through an intermediate triangle), the 
molecular orbital is significantly broadened. In the n = 1 
rectifier, the sulfur is separated from the ring by a single 
insulating -CH2- group. However, because of the sp 3 hy- 
bridization of carbon, the sulfur atom is out the plane of 
the ring. Then the separation between the sulfur and the 
ring is not large enough to prevent elongated p-orbitals 
of sulfur and carbon from direct overlapping. As a result, 
the electronic level is still significantly broadened. The 
situation changes radically for n = 2. Two insulating 
groups move the sulfur away from the ring by 4.3 A, so 
that direct overlap between the sulfur and the ring wave 
function becomes small. (From this point, insertion of 
new groups results in an exponential decrease of over- 
lap.) According to our calculations, for n = 2 the width 
is just ~ 10 meV, see Fig. ||. 

The numerical results illustrate the general formula for 
transmission probability T(E), Eq. (Q). For the present 
molecular diodes, one partial level width is much larger 
than the other, say Tl Fr. Then transmission at 
the resonance is T(E — Emo) — ^Tr/Fl- Thus, as 
decreases due to better insulation, the resonance gets 
narrow but higher. For the n = 2 diode, the resonance 
is narrow enough to be comparable with thermal, disor- 
der, and other types of broadening present in the sys- 
tem. Further increase of the short barrier is unnecessary, 
it will result only in reducing the molecular asymmetry 
and spoiling the rectification property. Our conclusion 
from this analysis is that two -CH2- groups on the short 
side is the optimal choice for rectification. One should 
add that this rule has been derived for thiol-terminated 
molecules that are self- assembled on gold or other no- 
ble metal. For other types of contact, this may change. 
For instance, when an LB film is deposited on a metal- 
lic surface, a larger (Van-der-Waals-like) gap may exist 
between the film and the metal. Such a gap will serve 
as an additional insulator, which may reduce the optimal 
length of insulating material on the molecule itself. 

(ii) Effect of the metal work function on the forward 
and reverse voltages. According to Eqs. (||) and A3), the 
forward and reverse voltages are directly proportional to 
the equilibrium energy difference A between the conduct- 
ing molecular level and the Fermi energy of the metal. We 
have studied this effect by varying the electrode Fermi en- 
ergies with respect to the molecular orbitals. The results 
are presented in Fig. |^ and in Table [ll]. The I-V char- 
acteristics all have similar shape but their forward and 
reverse voltages systematically increase with A, as ex- 
pected. Analysis of parameters from Table |TJ reveals that 
the ratios V F /A « 1.6, Vr/A « 2.6, and Vr/Vf ~ 1.6 
all remain approximately constant as functions of A for 
m = 2 and n = 6. Thus the ratio Vr/Vf is indepen- 
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FIG. 6. The effect of the metal work func- 

tion on the I-V characteristic of the molecular diode 
-S-(CH 2 )2-C 6 H4-(CH 2 )6-S-. 

dent of electrode material and is indeed a characteristic 
of the molecule asymmetry only. However, the relation 
Vr/Vf = f] w L r ight /L\ e [t is satisfied only approximately. 
For the (2,6) molecule, presented in Fig^, L r i g ht /^left = 3 
which is almost twice the voltage ratio. This is because 
the simple equations (||) and (||) do not take into ac- 
count the voltage drop on the central conductive unit. 
Clearly, some voltage always drops there, which effec- 
tively increases the lengths of both barriers and reduces 
the anisotropy. We discuss this issue in more detail in 
subsection (iii) below. 

The last two columns of Table || present data on cur- 
rent rectification. As a measure of the latter, we choose to 
compare currents at some positive and negative voltages 



W (eV) 


A (eV) 


V F (V) 


V R (V) 


V R /V F 


^op (V) 




4.0 


0.55 


0.90 


1.44 


1.60 


1.13 


56 


4.2 


0.75 


1.19 


1.93 


1.62 


1.41 


118 


4.5 


1.05 


1.66 


2.73 


1.64 


2.07 


173 


5.0 


1.55 


2.51 


4.01 


1.60 


3.02 


276 


5.2 


1.75 


2.77 


4.50 


1.63 


3.68 


305 


5.5 


2.05 


3.25 


5.32 


1.64 


4.25 


323 



TABLE II. Parameters of the I-V characteristics shown 
in Fig. | for m = 2 and n = 6. A = W — -Elumo is 
computed from the work function W and the molecular or- 
bital energy -Elumo = —3.45 eV. Vf and Vr are the cur- 
rent onset voltages under positive and negative bias, respec- 
tively. V op is some operating voltage between Vr and Vf- 
The last column represents typical current rectification ratios 
/+//_= I(+V op )/I(-V op ). 



FIG. 7. Dependence of the I-V characteristic on the asym- 
metry of the molecular diode HS-(CH2)2-C 6 H4-(CH 2 )„-SH 
for work function W = 4.0 eV. Notice how the Vf goes down 
while Vr goes up, systematically shifting the I-V curve to the 
left from the origin. 

with the same absolute value V op . Obviously, rectifica- 
tion is strongest if V op is chosen between Vf and Vr. In 
this case, in the positive direction there is already some 
appreciable current due to resonant tunneling through 
the molecular level. In the reverse direction, the tunnel- 
ing is still under barrier and therefore is exponentially 
reduced. This is the essence of our rectification mecha- 
nism. The values of the operating voltage V op shown in 
the table are roughly halfway between Vf and Vr. (Ap- 
proximately, V^p ~ 1.3Vp.) One can see that the current 
ratio is a steadily increasing function of A. The overall 
dependence is close to linear, with both end points being 
slightly off this general trend. 

An interesting problem is choosing an optimal value of 
A. As we have seen, large As result in better rectifica- 
tion of current and wider rectification voltage windows 
Vr — Vf- The latter fact is very important because all 
sorts of disorder in the system will tend to shrink the 
working window. Also, a particular circuit design may 
require the operating reverse voltage must be larger in 
absolute value than a certain voltage. In general, the 
absolute value Vr has to be as large as possible. Both 
these arguments favor a large A. On the other hand, 
large operating voltages imply higher power and higher 
electrostatic stress on the molecules, which are, of course, 
undesirable. Quantitative understanding of this trade off 
problem warrants a detailed investigation for a particular 
application. 

(iii) Effect of the length of the insulating barrier on the 
asymmetry of the I- V characteristic and on the overall 
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current magnitude. This effect is central to the present 
paper. By increasing the length of the longer barrier 
while keeping the shorter one fixed, one forces a larger 
portion of the voltage to drop on the longer barrier. This 
implies that the conducting threshold is reached at differ- 
ent bias voltages in the two opposite directions, as illus- 
trated in Fig. |^. This design principle provides us with 
the means to systematically increase the current rectifica- 
tion ratio and produce better molecular diodes. In Fig. |?|, 
we show the calculated FV characteristics of a series of 
molecular diodes (2, n) for work function W = 4.0 eV. 
The electrical parameters of the FV characteristics are 
summarized in Fable If] . The series begins with the sym- 
metric molecule n = 2, which produces a symmetric FV 
characteristic with Vr = Vf = F19 V. With increasing 
n, Vf monotonically decreases but Vr increases, result- 
ing in a systematic shift of the FV characteristic with 
respect to the origin. By n = 10, the reverse-to-forward 
voltage ratio reaches Vr/Vf = 2.18. The voltage ratio 
grows slower than oc n/2, suggested by Eqs. (§) and (||). 
It could be described by the relation 



Vr 
V f 



bright + 



L 



left 



(14) 



where aL is a correction that takes into account the volt- 
age drop on the central conducting unit. Using the data 
of Table [ffl|, one finds AL/L lcit = 2.5, 4.0, 4.0, and 4.9 
for n = 4, 6, 8, and 10, respectively. The correction does 
not remain constant as a function of n. Therefore it is 
difficult to assign to it a useful geometrical meaning. 



The penultimate column of Table III shows the cur- 
rent rectification ratios It grows approximately 
proportional to (n — 2), reaching a value of ~ 100 at 
n = 10. 

The important issue is the exponential growth of the 
diode resistance with n, see the last column of Table OF 
This result is, of course, expected because the transmis- 
sion as well as the current through the molecule is directly 
proportional to Tr, see Eqs. (Q) and (||) at » Tr. 
As the barrier length increases, the probability of under- 
barrier tunneling goes down exponentially fast, which is 
reflected in the numerical data. According to our cal- 
culations, addition of every extra pair of -(CH2)- groups 



n 


V F (V) 


Vr (V) 


V R /V F 


Vop (V) 






2 


1.19 


1.19 


1.00 


any 


1 


0.33 a 


4 


0.99 


1.28 


1.29 


1.13 


32 


2.40 


6 


0.90 


1.44 


1.60 


1.13 


56 


42.6 


8 


0.81 


1.61 


1.99 


1.13 


76 


836 


10 


0.80 


1.74 


2.18 


1.13 


104 


13,150 



TABLE III. Parameters of the I-V characteristics shown 
in Fig. ^. R is the resistance of one molecule at V = +V op - 
a resistance at V = 1.6V. 
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19 


27 


47 


85 
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(2,4) 


32 


52 


91 


138 


156 


154 


(2,6) 


56 


74 


173 


276 


374 


444 


(2,8) 


76 


92 


225 


N/C 


438 e 


N/C 


(2,10) 


104 


118 


290 


N/C 


556 e 


N/C 


Vop (V) 


1.13 


1.51 


2.08 


3.02 


3.50 


3.97 



TABLE IV. The rectification ratio I+/I- for six values of 
the wire work function and nine different molecules. The first 
column is the (m, n) numbers of the diodes. The top row is 
the work function of the electrodes in eV. The bottom row is 
the operating voltage V op in volts, at which the rectification 
ratio is reported. "N/C" stands for non-conclusive evidence 
due to numerical noise in calculating small currents under 
reverse bias. Best rectification is achieved for the (2,10) MR 
at W = 5.2 eV (in bold). The W value is close to the work 
function of gold. 
a at V = 1.32 V. 
b at V = 1.89 V. 
c at V = 2.46 V. 
d at V = 3.59 V. 
e at V = 3.02 V. 



increases the resistance by a factor of 17 ± 2. [Note that 
for (n = 2) — > (n = 4) the transition should not fol- 
low the same trend because at (n — 2) the molecule is 
symmetric and Tl = Tr.] Thus we have an unfavorable 
trade off problem where the rectification ratio improves 
linearly with n but the resistance and the time constant 
of the device worsens exponentially. This is an inherent 
feature of our mechanism and in some sense the price for 
its simplicity. When discussing the rectifier resistances, 
one should not forget that what matters is the total resis- 
tance of an electronic element. With a typical target size 
of (10 x 10) nm 2 , the number of molecules per junction 
is going to be of order of 1000. Then even n = 8 and 
10 diodes will have total resistances of 1-10 MSI, which 
amounts to a time constant ~ 10~ 10 s for a lpm long 
nanowire, and ~ 10 _6 s for a 1 cm long nanowire. 

The data from Table || suggests that rectification could 
be improved further by increasing the initial energy dif- 
ference A. We computed I-V characteristics for every 
value of the work function listed in Table || and for all 
MRs (m, n) with m = 1, 2 and n = 2, 4, 6, 8, 10. [In addi- 
tion, everything was calculated for the two values of the 
electrode matrix element W ssa , bringing the total num- 
ber of studied FVs to 120J All the curves have shapes 
similar to those in Figs, g and 0, therefore we do not 
show them explicitly. Table |y| summarizes the results 
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(c) 




FIG. 8. (a) A gap in the density of states of the electrodes, 
(b) Under forward bias, when the right Fermi level aligns 
with the LUMO, there are still states available in the left 
electrode. The current flows, (c) Under reverse bias, by the 
time the left Fermi level aligns with the LUMO, there are no 
states available in the right electrode. The elastic current is 
blocked. 

on rectification ratios for a variety of combinations. As 
expected, the rectification ratio grows with both W and 
n. The maximum ratio of above 500 was observed for 
n = 10 and W = 5.2 eV (A = 1.75 eV). 

(iv) Role of energy gaps in the electrodes' density of 
states. This effect is not generic and potentially applies to 
semiconductor electrodes only. Nevertheless, it offers an 
interesting possibility to enhance the rectification prop- 
erty of the diode by inhibiting current flow in the reverse 
direction due to an energy gap in the density of states. 
The basic idea is illustrated in Fig. ||. At some voltage 
V g , the bottom of the gap aligns with the conducting 
molecular orbital. At higher V > V g resonant tunneling 
becomes impossible due to the lack of available final elec- 
tronic states. Accordingly, the current must drop to zero 
at V = V g . ^ V g lies in between Vp and Vr, then there 
will be some substantial current in the forward direction 
in the interval Vp <V <V g , Fig. |](b), while no current 
at all in the reverse direction, Fig. g(c). In the latter 
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-6 -4 -2 2 4 6 

Bias voltage (V) 

FIG. 9. I-V characteristic of the (2,8) MR for the electrode 
matrix element W ss <j = 0.909 eV. A gap in the electrode's 
density of states cuts off current at large voltages. Apart 
from current rectification, this leads to a significant negative 
differential resistance effect. The latter will be reduced if 
inelastic processes are taken into account. 

case, the gap of the right electrode is reached before the 
left Fermi energy aligns with the conducting orbital. 

The above argument does not take into account inelas- 
tic processes. If included, they would result in a non-zero 
current even in the situation of Fig. ||(c). In this case, 
the electrons can tunnel under the right barrier with irra- 
diation of phonons, arriving at the right electrode with a 
reduced energy for which there are available states at the 
top of the valence band. In this paper we assume that the 
clastic tunneling is the dominant transport mechanism 
and such processes can be neglected. We have modeled 
this effect by using the smaller electrode matrix element 
W ss <7 = 0.909 eV. Typical I-V characteristics are shown 
in Fig.^j and their parameters are summarized in Table 0. 
With increasing W (and A), the forward and reverse volt- 
ages Vp and Vr increase in accordance with Eqs. (||) and 
([|). At the same time, V g decreases because the bottom 
of the energy gap moves closer to the conducting molec- 
ular level. As a result, the window of resonant tunneling 
is squeezed from both sides and at some critical W van- 
ishes altogether. This effect can be seen in Fig|. In the 
reverse direction, the peak disappears between W = 4.2 
eV and W = 4.5 eV. In the forward direction, the same 
happens between W = 5.2 eV and W = 5.5 eV. (Of 
course, resonant tunneling survives longer in the forward 
direction because Vp > Vr.) 

Notice that the TVs of Fig. ^ display large negative dif- 
ferential resistance under positive bias. In this respect, 
our molecular diodes behave similarly to the conventional 
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semiconducting tunneling diodes. However, under nega- 
tive bias, the behavior is quite different. While the tun- 
neling diodes conduct very well under reverse bias, our 
MR conduct very little. 

V. CONCLUSIONS 

The major purpose of this paper has been to predict 
several trends in the I-V characteristics of MR when they 
are dominated by resonant tunneling through a certain 
molecular orbital localized on the conjugated part of a 
spatially asymmetric molecule. To achieve a large recti- 
fication ratio, the conjugated part of the molecule must 
be connected to electrodes by insulating molecular groups 
of different length. By varying the ratio between the bar- 
rier lengths, one can achieve a rectification ratio of several 
hundred while keeping the current through the molecule 
at measurable levels. 

As far as the described mechanism of molecular rec- 
tification is concerned, the observed trends could be di- 
vided into "spatial" and "energetic" domains. The spa- 
tial trends refer to the changes in the I-V characteristic 
that follow from changing the lengths of the two insulat- 
ing barriers of the MR. The length of the shorter barrier 
controls the width the molecular resonance. We have 
found the optimal length to be two -(CH2)- groups. One 
such group does not provide sufficient insulation of the 
conducting unit from the metal. As a result, the reso- 
nance is too broad and current asymmetry is not pro- 
nounced. More than two -(CH 2 )- groups make the res- 
onance narrower than temperature and disorder induced 
widths, which is not useful. Increasing the length of 
the longer barrier improves the current rectification ratio 
but increases the total resistance of the molecular diode. 
These trends can be observed experimentally by study- 
ing several molecules with different lengths of the longer 
chain but identical otherwise. The energetic trends re- 
fer to the dependence of the rectification property on the 



W (eV) 


A (eV) 


V F (V) 


V R (V) 


V s (V) 


V op (V) 


I+/I- 


4.0 


0.55 


0.84 


1.62 


N/C 


1.23 


70 


4.2 


0.75 


1.17 


2.16 


N/C 


1.80 


151 


4.5 


1.05 


1.58 


N/C 


5.20 


2.17 


282 


5.0 


1.55 


2.33 


N/C 


3.68 


2.83 


1308 


5.2 


1.75 


2.60 


N/C 


3.10 


2.83 


1064 


5.5 


2.05 


N/C 


N/C 


N/C 


2.83 


3 



TABLE V. Parameters of the I-V characteristics of the 
(2,8) molecular rectifier for W ssa — 0.909 eV. (I-Vs are shown 
in Fig. ^) Notice how the rectification ratio exceeds 1000 
when the resonant conduction disappears in the reverse di- 
rection, but then drops to just 3 when it disappears in the 
forward direction too. 
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FIG. 10. Other possible molecular diodes. 

work function of the electrode material and on the elec- 
tron affinity of the molecule. The two quantities define 
the parameter A, which is the most important energy 
scale in the problem. In general, larger A leads to bet- 
ter rectification but at the same time to larger operating 
voltages. These trends can also be checked experimen- 
tally by measuring, for instance, the same molecule on 
gold and silver electrodes, or by measuring the molecules 
with different conducting units (say, benzene and naph- 
thalene) on the same set of electrodes. 

Most of our results are rather insensitive to changing 
the end ("anchor") groups. The role of the latter is to 
provide better connection to the electrodes. Thus the end 
groups could be adjusted to the particular experimental 
setup without significantly altering the electrical proper- 
ties of the device. They could even be absent altogether 
if the measurement technique does not require covalent 
bonding between the molecules and the electrodes. The 
composition of the insulating barriers may be changed 
too. For instance, instead of alkane chains one could use a 
combination of the alkane segments and saturated cyclic 
hydrocarbons, see Fig. |l^. An advantage of this design 
is that the molecules have more or less the same cross 
section along their entire length. Therefore the packing 
of the film will be much better. On the other hand, such 
MR could be harder to synthesize. 

These conclusions are certainly limited to the partic- 
ular set of approximations used in the present study. 
For instance, we have ignored inelastic processes, which 
should become progressively more important for larger 
molecules and at larger bias voltages. Inelastic scat- 
tering generally increases the current through the sys- 
tem, since the carrier has more final states to tunnel 
into when it can lose energy by exciting a vibronic ex- 
citation. In the present resonant case, when V exceeds 
the threshold, either Vp or Vr, there will be some elec- 
trons with energy higher than the energy of the LUMO 
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resonance. Such electrons can still tunnel through the 
junction effectively by first loosing their energy and then 
tunneling resonantly. Thus the overall current will con- 
tinue to grow after the initial sharp increase, in con- 
trast to the elastic-only case where the current stays ap- 
proximately constant. Additionally, a finite temperature 
should smoothcn the sharp features of the I-V charac- 
teristic, thereby reducing the rectification ratio. The in- 
clusion of inelastic processes would probably reduce the 
estimated rectification ratio. 

In conclusion, we have described a simple principle for 
molecular rectification. It requires only one conducting 
molecular level located asymmetrically with respect to 
the electrodes. Theoretical calculations have shown that 
current rectification ratios of several hundred are achiev- 
able with this mechanism. Higher ratios are likely to be 
accompanied by a significant reduction in conductance, 
which will make such diodes less practical. We have 
made several predictions of how the rectification prop- 
erties should change between different molecules and dif- 
ferent electrode materials (work functions). These trends 
are experimentally verifiable on families of structurally 
similar molecules. Such trends are signatures of resonant 
tunneling through molecules and could serve as proof of 
experimental observation of the latter. 
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